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ABSTRACT p_ 

This Report summarizes the analyses of the interaction of the Sur- 
veyor radar altimeter doppler velocity sensor (RADVS) system and 
vernier thrust chamber plumes. The investigation was occasioned by a 
study task initiated by JPL. 

In keeping with guidelines of the study task outlines, this analysis 
was specifically conducted to provide conservative estimates of the 
vernier plume, RADVS-system interactions effects. In all instances, 
operating conditions were assumed which, while consistent with Sur- 
veyor system specifications, would exaggerate the RADVS-plume inter- 
action. For example, the vernier propellant specifications permit a 
maximum total impurity concentration of 500 ppm. Actually, the pro- 
pellant purity level is known, from test, to be at least an order of 
magnitude less than the specifications. However, for the purpose of 
this investigation, the propellant was assumed to have 500 ppm of 
sodium as an impurity. As with all other assumptions made during the 
course of analysis, this assumption permitted an upper bound to be 
placed on the interaction of the vernier plumes and the RADVS system. 
As a consequence, no attempt was made to provide a detailed assess- 
ment of the most probable level of RADVS-plume interaction to be ex- 
pected under nominal Surveyor operating conditions. 

The maximum calculated two-way attenuation and reflection level 
are both well below the levels at which the performance of the RADVS 
system would be significantly influenced. However, the nonstatic two- 
way phase-shift was found to be in excess of the level above which 
the operation of the RADVS system would be affected. 
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I. INTRODUCTION 


The Surveyor spacecraft has been expressly designed 
to permit the controlled deposition of a scientific and/or 
engineering payload on the lunar surface. This feat is 
to be accomplished by means of a controlled retro-rocket 
descent, from an Earth transfer trajectory, to within 
closely prescribed distances from the Moon. 

Basically, the retro sequence conducted in the vicinity 
of the Moon can be divided into three distinct phases 
(Ref. 1). During the first phase, a radio altimeter located 
in the nozzle of the main solid retro-rocket provides an 
altitude ( — 100 km) marker signal which, when processed 
by the control circuitry of the spacecraft, initiates the 
main retro-rocket sequence. During retro bum, the three 
bipropellant vernier thrust chambers, used also for mid- 
course trajectory corrections, provide attitude control. 
Upon completion of the main retro-rocket burn, the 
spacecraft is within —12 km of the lunar surface, and 
its velocity has been reduced from —2700 m/sec to 
— 100 m/sec. Twelve seconds after completion of the 
main retro burn, the solid propellant engine is ejected. 
At this point the spacecraft is —8 km from the surface 
and descending at a relative velocity of —110 m/sec. The 
initial phase of the landing maneuver is then complete. 

The spacecraft is then at a sufficiently low altitude, 
relative to the surface of the Moon, that meaningful 
signals can be received by the doppler antennas and the 


radio altimeter. The data from these antennas, which 
form a part of the radar altimeter and doppler velocity 
sensor (RAD VS) system, are used to control the rate 
of descent and the attitude of the spacecraft by means of 
throttling of three vernier engines. Pitch and yaw are 
controlled by differential throttling of the verniers, while 
roll control (about the vertical, z-axis, of the spacecraft) 
is maintained by gimbaling (hinging) one of the vernier 
thrust chambers (vernier thrust chamber No. 1). The 
vernier engines operate continuously throughout the 
descent phase until an altitude of — 4 m from the lunar 
surface is reached. At this point, if the flight control 
electronics have functioned as planned, the remaining 
velocity components, relative to the surface, will be less 
than —2 m/sec and the vernier engine thrusting is 
terminated. 

During the third, and last, phase of the descent the 
spacecraft free-falls to the surface, the impact being 
absorbed by the structure of the bus. 

Early in the Surveyor program, a study was conducted 
to determine whether, during the second phase of the 
landing sequence, the RADVS system would suffer a 
degradation in performance because of the free electrons 
present in the vernier exhaust plumes, through which the 
doppler velocity and radio-altimeter beams propagate. 
The analysis conducted by P. Molmud (Ref. 2) of STL 
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indicated that there would be no appreciable effect of the 
exhaust plumes on the RADVS system, for a similar space- 
craft system. 

This study is a reevaluation of the degree of interaction 
to be anticipated between the RADVS and the vernier 


exhaust plumes. More recent data than were available to 
Molmud concerning the physical characteristics of the 
spacecraft have been used. In addition, a more careful 
evaluation of the anticipated back-scatter and the effects 
of throttling and gimbaling of the vernier engines has 
been included. 


II. SPACECRAFT GEOMETRY 


As seen in the plan view of Fig. 1, the two RADVS an- 
tenas and the three vernier thrust chambers are arranged 


on a circle that is roughly 1.9 m in diam. The exit plane of 
the vernier engines exhaust nozzles describe a plane 



-y 

Fig. 1. Plan view Surveyor 
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I 


which, to within narrow limits, passes through the center 
of the RADVS antenna apertures. The plan view depicts 
the directions in which the various doppler and radio 
altimeter beams propagate. A clearer picture of the 
geometry of these beams with reference to the spacecraft 
is presented in Fig. 2. All of the doppler beams lie, ap- 
proximately, in a cone (half angle of 25 deg) with respect 
to the vertical spacecraft axis (z axis). 

Consideration of these diagrams leads one to conclude 
that the most severe beam-plume interaction would occur 
between doppler beam No. 2 and the plume from vernier 
thrust chamber No. 1. This conclusion is based on the fact 
that the half of the antenna generating beam No. 2 and 


the vernier thrust chamber No. 1 are in closer proximity 
than any other combination of antenna and thrust 
chamber and, moreover, vernier thrust chamber No. 1 
is gimbaled to provide spacecraft roll control. This 
gimbaling would, in effect, permit the plume to intersect 
more or less of the beam— depending on hinge position— 
therefore, changing the electron density environment 
through which the beam propagates in a variable manner 
over and above that which can be expected solely from 
the effects of throttling. Therefore, since the interaction 
of vernier thrust chamber No. 1 plume and the doppler 
beam No. 2 is the most severe, both statically (no 
throttling, no gimbaling) and dynamically (throttling and 
gimbaling of the thrust chamber), the study will be re- 
stricted to this case to the exclusion of all others. 


-z 



Fig. 2. RADVS beam orientation 
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III. VERNIER THRUST CHAMBER NO. 1 PLUME, 
DOPPLER BEAM NO. 2 GEOMETRY 


Figure 3 presents a diagram of the geometric relation- 
ships that exist between the center line of the plume 
generated by vernier thrust chamber No. 1 and the beam 
pattern produced by doppler antenna No. 2. From scale 
drawings of the spacecraft, the inter-axis spacing of 
vernier thrust chamber No. 1 and doppler antenna No. 2, 
a = 42 cm. The vernier thrust chamber gimbals in a plane 
that intersects the plane formed by lines parallel to the 
z axis (drawn through the thrust chamber and antenna), 
at an angle /3 — 12.7 deg. The gimbal angle y can vary 
over the range of ±5Y> deg. For this analysis, only two 


cases will be considered; the ungimbaled case where the 
thrust axis is vertical and y = 0, and the case where 
thrust chamber No. 1 is gimbaled the maximum amount 
( 5Vj deg) in the direction toward doppler beam No. 2. 
We shall denote distances along the plume axis (as 
measured from the thrust chamber exit plane) as $ and 
distances measured from the antenna aperture plane 
along the beam axis as rj. The first requirement is to deter- 
mine the distance of closest approach (R) of a given beam 
ray to an arbitrary point | on the plume axis by consider- 
ing the geometry of the antenna beam. 


THRUST 



Fig. 3. Vernier thrust chamber No. 1 , doppler beam 
No. 2 geometry 
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IV. ANTENNA BEAM GEOMETRY 


Two elliptical antennas, which are split so as to gen- 
erate two separate and distinct beams for a total of four, 
comprise the antenna subsystem of the RADVS. The 
individual beams emanating from one half the RADVS 
antennas can be considered as being produced by an 
antenna having dimensions of ~~25 X 38 cm. Using the 
average of these dimensions to characterize the antenna 
aperture dimensions, i.e., D ~ 32 cm, the calculated 
Rayleigh range l = D-/2A. 0 , for a frequency of 13.3 GH* 
(A„ = 2.25 cm), is ^2.3 m. This implies that the beam will 
be approximately collimated for ^ ^ 2.3 m and then 
diverge at an angle a, where tan a — A 0 /D, or a ~ 4 deg 
for 77^2.3 m. Two beam rays are of particular interest, 


the center ray representative of maximum beam power 
and the — 3-db (half power) ray closest to the plume. 
The —3-db ray will be assumed to be emanating from the 
edge of the antenna, i.e., at a distance of D/2 from 
the central (maximum power) ray. The —3-db ray is, 
in the Rayleigh near field region 77 taken as parallel 
to the center ray and to lie on the surface of a cone of 
half-angle a ~ 4 deg for *7 ^ Z, i.e., outside the Rayleigh 
near field zone. As will be subsequently demonstrated, a 
preponderant portion of the beam-plume interaction will 
occur in the near field region of the antenna, the fact that 
the assumed —3-db ray diverges for 77 is of relatively 
little importance, but is included for completeness. 


V. PLUME GEOMETRY 


A . Gas Density 

It will be assumed that there is no interaction of beam 
No. 2 with other than electrons present in the plume 
from vernier thrust chamber No. 1. This assumption is 
deemed justified on the basis of the thrust -chamber spac- 
ing and the overall geometry of the plume. The plume 
geometry assumed for this program is that calculated 1 to 
exist for a conical (15-deg half-cone angle) nozzle ex- 
hausting into a vacuum, a gas with a ratio of specific heats 
C p /Cr — 1.313. This ratio of specific heat is equivalent to 
what one would expect for the propellants used, i.e., 
nitrogen tetroxide (with 10% nitric oxide) and mono- 
methylhydrazine monohydrate. While the actual vernier 
thrust chambers do have contoured nozzles, the error 
introduced into the calculation of the electron density 
profiles in the plume by the assumption of a 15-deg half- 
angle cone nozzle is considered to be insignificant. 

In Fig. 4 are shown the Mach number profiles for the 
assumed conical nozzle, plotted in terms of dimensionless 
radial and longitudinal coordinates, i.e., coordinates 
normalized in terms of the thrust chamber exhaust plane 


'McDermitt, C. Eugene, to Briglio, Jr., A., “Impingement Heating 
Rate Estimates of the Current R\1D Exhaust Plume in the Vicinity 
of the Crushable Blocks and Altitude Marking Radar,” Private 
Communication, Jet Propulsion Laboratory, Pasadena, Calif., 16 
Sept. 1965. 


radius (R e = 6.5 cm). Also shown in the figure are the 
beam rays calculated for the four cases of interest, i.e., 
two gimbal positions of vernier thrust chamber No. 1 
(y — 0 and 5 l /*>-deg) and for each of the two beam rays 
of interest (central and —3-db ray closest to the plume). 
It should be noted that even though the rays have been 
assumed straight (this will be substantiated later) they 
appear curved in the figure. This is because the beam 
rays and the center line of the plume are askew, the 
radial distance R between the two being a non-linear 
function of the plume coordinate £. 

The isentropic perfect gas relation 

( k — 1 

1 + — Mi xit 

1. i 

1 4 £ local 

can now be used to evaluate the local gas density p as a 
function of the known Mach number M contours and the 
nozzle exit conditions. The following nozzle exit condi- 
tions are those anticipated to exist for the maximum thrust 
(47 kgf) condition 

10 -6 g 

Pexit — 4.04 X —~T > 



5 



JPL TECHNICAL REPORT NO. 32-927 



Fig. 4. Plume, RADVS beam geometry 


( N a — 1.07 X 10 17 cm 3 for (}]7exit — 22.77), where N a is the 
heavy particle density and 07 ? is the average molecular 
weight, 

M„ xit = 5.20, k = ^ = 1.313 
T exit = 840°K 

The nominal expansion area ratio for this thrust chamber 
is 86:1. 


B. Electron Density N t 

To evaluate properly the electron density level present 
in the plume, it is necessary to consider the ionization 
processes associated with combustion. This, in turn, re- 
quires a knowledge of the level of low ionization level 


impurities (alkali metals, etc.) present in the combustion 
gases. There are two sources of possible contamination, 
the propellants per se, and the material ablated from the 
chamber, throat, and nozzle liners. The JPL propellant 
specifications GMO-50427-MAT-A and GMO-50428- 
MAT-A require that there should be less than 500 ppm 
total inpurities present in any batch of propellant speci- 
fied for use. No information is currently available con- 
cerning the level of contamination caused by impurities 
present in the propulsion flow system which could be 
absorbed by the propellants during extended periods of 
storage. This source of contamination is thought to be 
small, i.e., less than 10 ppm. Actual analysis of samples 
of propellant indicate that extraneous impurities are less 
than 10 ppm, i.e., well below specification. The Rokide Z 
and silicon carbide coatings on the internal surfaces of 
the thrust chamber will erode slightly during operation. 
At most, it is estimated that — 6 g of material will be 
ablated. The total propellant flowing through the system 
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is '-'84 kg. Therefore, ablation could contribute, on the 
average, a maximum of 60 ppm of contamination. 

To cover all possible contingencies, it will be assumed 
that there are 500 ppm of contamination present in the 
exhaust gases. Moreover, to exaggerate the effect of this 
contamination, it will be taken to be comprised solely of 
sodium. These assumptions are quite conservative, insofar 
as the alkali metal content in the propellant is most prob- 
ably on the order of 10 ppm or less. Consequently, the 
electron density, since it arises almost exclusively from 
ionization of the alkali metals present in the chamber 
gases, will be overestimated by at least an order of 
magnitude. 

The chamber conditions used in the evaluation of the 
ionization level were: P c — 17 atm, N a — 4.2 X 10 ,9 em~\ 
T c = 2950° K, gas composition — combustion products 
having an ionization potential of '—11 ev plus 500 ppm 
of sodium (ionization potential 5.1 ev). 

Application of the Saha equation for this two-species 
gas showed that only a negligibly small portion of the 
ionization level was due to the combustion products gas, 
the indicated ratio of electron to atom density 

= c~2X 10* 

N a 

arising from ionization of the sodium. 


approximately 200 ppm (there were no tests at the 500 
ppm level assumed in this study) and chamber conditions 
comparable to those encountered in the Surveyor vernier 
system, the electron density was, to within experimental 
error (a factor of ^2), equal to that calculated for the 
exhaust plane of the vernier thrust chambers using the 
above assumptions. 

It should be noted that the electron density is de- 
pendent upon the square root of the sodium concentration 
and that the thrust chambers tested had a somewhat 
smaller expansion area ratio than those used on Surveyor. 
Nevertheless, the agreement between the test results and 
the above calculations is exceptional and validates the 
assumptions used. 

The assumption of a fixed ionization level, frozen at 
the averaged chamber-throat conditions, exaggerates the 
electron density level and is, therefore, conservative. This 
assumption permits one to relate the electron density level 
N e very simply to the heavy particle density N a , viz. 
N e — c N a . For the ensuing calculations, it is convenient 
to normalize the electron density to the electron density 
level (N p ), which would just cause cutoff of the doppler 
beam signal. The doppler beam is generated at a fre- 
quency of 13.3 GHz. Therefore, equating the beam 
frequency to the plasma frequency, 

*> =2 7T X 13.3 X 10* = a* = 


£ofn e ) 


The ionization level calculated on the basis of chamber 
conditions was ^3X 10“ and on the basis of throat con- 
ditions, —0.7 X I0~ r \ It is anticipated that equilibrium 
conditions would prevail up to the throat and that the 
flow would be frozen downstream of this station. Accord- 
ingly, the ionization level was taken to be —2 X 10 “. 
For comparison, if the flow were assumed to be in 
equilibrium through the nozzle and to be frozen in the 
plume the ionization level would be completely negligible, 
i.e., less than 10~ 10 . 

A series of tests have been conducted by W. Balwanz 
of the Naval Research Laboratories 2 (NRL) to determine 
exit plane electron densities of 23 kgf thrust chambers on 
alkali metal seeded Surveyor vernier thrust chamber pro- 
pellants. The X-band focused microwave probe used in 
these tests indicated that, with sodium concentrations of 

3 The results of these tests were communicated by H. Schuman of 
Hughes Aircraft Co. H. Schuman assisted Dr. Balwanz in con- 
ducting these tests at NRL. The results, as yet, have not been 
formally documented. 


where e and m e are the electronic charge and mass, re- 
spectively, and e Q is the permittivity of free space, we find 
that: N p = 2.2 X 10 1 - cmr'\ and where N p is the cutoff 
electron density. 

Figure 5 is a plot of the normalized electron density 
N e /N„ vs the distance along the antenna-ray axis ? j normal- 
ized to the freespace wave length A 0 — 2.25 cm. The four 
cases of concern are shown, i.e., the two rays and two 
gimbal angles of interest. As indicated previously, the 
bulk of the beam-plume intersection will occur within 
the near field of the antenna (inside the Rayleigh range, l) 
by virtue of the fact that most of the electron density 
encountered along the ray path is concentrated well 
within this range. 

C. Electron Collision Frequency, v c 

The attenuation of an electromagnetic wave propagat- 
ing through a plasma (Ref. 3) is dependent on the degree 
of damping the individual electrons undergo when forced 
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l.5[ 


/ \ DOPPLER BEAM 2, VERNIER THRUST 


CHAMBER I PLUME 


CENTER RAY 


\ 3 db RAY CLOSEST TO PLUME 

\ N e /N a = c CONSTANT = l.75xlO -6 BASED “| 
\ ON CHAMBER CONDITIONS £=2950 "K, 
P c -\ 7 atm, 500 ppm SODIUM 

RAYLEI6H RANGE; £> = 32 cm X 0 =2.25cm| 
f=l3.3GHz, A£ = 2.2x|0 12 cm' 3 



y o 5.5 deg, £ = 12.7 deg THRUST 
CHAMBER I GIMBALED IN 
DIRECTION OF BEAM 2 


200 


NORMALIZED POSITION ALONG BEAM 2 


Fig. 5. N e /N p vs r)/\ Q 


into oscillation by the electric field of the wave. At the 
low ionization levels encountered in the plume these 
collisions will, almost exclusively, be with neutral atoms. 
The averaged electron momentum exchange collision 
rate is: 




depends upon the square root of temperature, this dis- 
tinction could lead to a factor of —2 increase in v e . 
Perusal of published elastic electron-atom momentum 
exchange cross section data for gases, representative of 
the combustion products and the alkali metals, permitted 
a reasonable estimate of the maximum expected cross 
section Q ea to be made. This value is Q ea = 3 X 10 ir> cm 2 . 
Now for T e ~ 2950°K, v e ~ 3 X 10 7 cm/sec; 


v r — 9 X 10 N N a • 


As we have in Fig. 5 a plot of the normalized electron 
density along the various rays of interest, it is expeditious 
to relate collision frequency to N e /N p and additionally, to 
normalize the frequency by the beam frequency w. 



so: 


Vy 

O) 


= 9 X 10 8 




C to 



where v e is the electron mean thermal velocity and Q ea 
the average electron atom cross section for elastic mo- 
mentum exchange. As we wish to be conservative, i.e., 
determine the maximum possible collision frequency and 
hence, attenuation, the chamber temperature T e will be 
used to calculate v e as opposed to the use of the exit tem- 
perature, which is perhaps more realistic. As the velocity 


As seen from Fig. 5, 

< 2 X 10 1 — < 3 X 10-* 

Np <*> 

therefore, the plasma only very weakly attenuates the 
beam. 
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VI. BEAM-PLUME INTERACTION 


A. Relative Propagation Constant K* 

The propagation of a plane, plane-polarized monochromic electromagnetic wave 
in an isotopic inhomogeneous plasma (inhomogeneous in the direction of propaga- 
tion > 7 ), can be described in terms of its electric (£) field by the wave equation: 

d-E , ~ 

w +kaME =° 

where E = Ee tta \ and k* ( 77 ) is the complex propagation constant. It is convenient 
to normalize k* by the propagation constant for free space (k 0 — 2 ? t/\ 0 ). 



It can be shown that K* (Ref. 4) is related to the macroscopic properties of the 
plasma, i.e., N e and v, as follows: 

K ‘ = m{~ < l - P) + [<* - p > ! + f ‘(t)7T 


where 


N e 



For the present case N e /N p and v r /<o < < 1, which permits us to write: 



where of course 


We = NeM 
N p N p ' 

Note that we could write the imaginary part of the relative propagation constant 
K, as 0.7 ( N e /N p ) 2 where we have assumed T c and Q af (and c) as being constant. 
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B. Applicability of the Phase Integral ( WKBJ) Method 

Where gradients in the propagation constant are small, it is possible to use the 
phase integral approach, i.e., geometrical optics, to describe the propagation of a 
wave through an inhomogeneous media. In general, this technique requires that 
the electromagnetic wave be plane, i.e., the direction and amplitude are the same 
over distances, which are large in comparison to a wave length, and that the 
direction of propagation is normal to the wave surface. In essence, geometrical 
optics correspond to the limiting case of \ 0 —> 0. The criteria that needs to be 
satisfied for the application of this approximation is: 


VK* _ % r 


k* ' v<< x„ 


where ?/ is the unit vector in the direction of propagation. 


Substitution of the expressions for K * into the above relationship results in: 


4t) .H® 


<< 2t> 


This criteria is easily satisfied even for the worst case of the — 3-db ray, y = 5V> deg 
situation at the point of maximum slope 10, cf Fig. 5). It should be pointed 

out that the near field region of an antenna is generally characterized by quite 
complex (non-planar) phase fronts. In addition, the wave is propagated through 
a three-dimensional electron density profile, and we have tacitly assumed that the 
plasma was uniform, transverse to the direction of propagation. While these two 
factors would tend to vitiate the use of geometrical optics, their precise influence 
is thought to be small. That is, uncertainties associated with the use of this tech- 
nique are commensurate with those introduced by the previous assumptions. 


The geometrical optical or phase integral technique permits one to describe the 
propagation of a wave in an inhomogeneous medium in terms of local conditions 
along the ray path 


E 



K* dr] 




C. Phase Shift 

The phase <j> p of the wave propagating through the plasma is dependent upon 
the real part of the propagation constant, i.e.. 
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The phase change A <j> is then the difference between the phase in plasma and 
that calculated for free space conditions, i.e., <j> a = 27rA(iy/A rt ) 

Therefore, 



for one-way propagation. 

The two-way phase shift due to the presence of the plume and exclusive of phase 
shifts caused upon reflection from the surface is: 




rad 


The normalized electron density profiles of Fig. 5 were developed on the basis 
of maximum thrust (47 kgf) operation of the thrust chamber. To evaluate the 
effects of throttling upon the phase shift, it wiii be assumed that the combustion 
efficiency and flow characteristics of the nozzle are independent of thrust level. 
With these assumptions, then N e /N p would be directly proportional to the thrust 
level. Table 1 summarizes the calculated (by numerical integration of the N e /N p 
profiles) two-way phase shifts to be expected for the two ray paths, two gimbal 
positions and two thrust (T) levels (maximum — 45 kgf, minimum = 13.6 kgf). 


Table 1 . Two-way phase shifts 


Condition 

y = 0 deg 

y = 5Vj deg 

This* 

T ID i II 

Timx 

T,„ j i, 

Central ray, deg 

-7.3 

-2.1 

-16.3 

-4.2 

— 3-db ray closest 
to plume, deg 

-9.9 

-2.8 

1 

-24.1 

-7.1 


The doppler velocity sensing system would be affected by variable phase shifts. 
Therefore, it is meaningful to present the data of Table 1 in a somewhat different 
format, which would emphasize the effects that gimbaling and throttling of the 
vernier thrust chamber would have on the two-way phase shift. To do this, we 
take as a nominal condition, y = 0, and an average thrust level of 30.4 kgf. This 
corresponds to a phase shift of —4.7 deg for the central ray and —6.3 deg for 
the — 3-db ray closest to the plume. These calculations are summarized in Table 2. 


Table 2. Change in two-way phase shift from normal 


Condition 

y = 0 deg 

y = 5 Vi deg 

Tmsix 

Tm in 

Timx 

T„,in 

Central roy, deg 

-2.6 

+ 2.6 

— 11.6 

+ 0.5 

— 3-db ray closest 
to plume, deg 

-3.6 

+ 3.5 

-17.8 

-0.8 
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The preceding calculations were based on estimates of the maximum possible 
plume electron density levels. Analysis of the RADVS system suggests nonstatic 
phase-shifts less than one deg would not significantly influence overall system 
performance. The actual electron density levels and the phase shifts are probably 
lower than those shown in Table 2. However, the results indicate that the plume 
could interfere with the doppler velocity sensor system. 


D. Attenuation 

The attenuation of the beam is dependent upon the imaginary part of the rela- 
tive propagation constant, i.e., 


“p 1 ■ K ' (’)] d" = C. “ p [ " ’ (^rXir )] d (i) 


or the attenuation of the electric (E) field strength for one-way propagation of 
the wave is: 



exp 


[-«■ (tTKi) 


Since the signal strength is proportional to E~ the overall attenuation of the signal 
because of two-way propagation through the plume, neglecting the losses asso- 
ciated with the reflection process is: 


0* p [“ 5 Kx) Mt) 


In terms of db of attenuation, i.e., 


db attenuation — 101og to y* exp 5.6?r J d 

The results of the numerical integration of the above expression for the two rays 
and the two gimbal positions of interest are shown in Table 3 for the case of 
maximum thrust. 



Table 3. Two-way signal attenuation 


Condition 



y = 0 deg 

y = 5 Va deg 

Central ray, db 

— 3-db ray closest 
to plume, db 

“3.2 X 10 4 
-4.9 X 10 4 

-1.7 X 10 3 
-4.0 X 10 :t 


No attempt was made to determine the attenuation associated with the less 
severe minimum thrust case. In all cases the indicated attenuation is well below 
0.1 db, a value at which attenuation effects begin to noticeably affect the per- 
formance of the RADVS system. 
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E. Reflections from the Plume 

Application of the phase integral method to the problem of reflections from 
arbitrary electron density profiles shows that in the first approximation there 
is no reflection. However, the following expression for the reflection coefficient 
X R (i.e., the ratio reflected to incident E fields) is obtained (Ref. 3) using the 
second order terms of the series 



To facilitate computation, the small amount of damping represented by the 
imaginary part of K *, i.e., K h will be neglected. This permits us to write: 





having set q Q / A„ equal to zero. Because of the involved nature of this expression, 
only one case was numerically evaluated: the y = 5 1/2 deg, — 3-db ray case. The 
quantity of ultimate interest is the reflection coefficient for the signal, defined 
as tRRR*, where /JP is the complex conjugate of fR. The results of the numerical 
evaluation of this quantity indicate that the reflected energy would be a factor 
of —2.2 X 10 11 of the transmitted energy. In terms of db then, the reflected energy 
is — —106 db below the transmitted signal. As the operation of the RADVS 
system would not be affected by the signals reflected from the plume which are 
— 90 db, or more, below the level of the transmitted signal, the presence of the 
plume would not impair the proper operation of the system. 


F. Bending of the Beam 

The radius of curvature (R c ) of the beam ray can be shown (Ref. 4) to be 
represented by the expression: 



However, implicit in our use of a basically one dimensional model of the plume, as 
discussed above, is the assumption that the plume is uniform in a direction normal 
to the ray (direction of propagation). It can be rigorously demonstrated (Ref. 4) 
that no matter how involved the variation of the propagation constant is in a region 
bounded by two planes q = 0 , — A), the change in ray direction will be 

determined by the propagation constants to the left of q = 0 and to the right of 
q ~ A, provided N e < N p . In our case, the region close to the antenna (77 — 0) and 
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far from the plume (77 ^ A) are regions characterized by the vacuum propagation 
constant, k 0 . Therefore, use of the simple (one dimensional) model leads to the 
conclusion that there will be no bending of the ray. However, since it was not 
possible in the time allocated to this study to calculate the effects of the three 
dimensionality of the electron density profile, it was necessary to try and make 
an estimate of the maximum possible curvature the beam could undergo. This 
was done by assuming that the maximum gradient in electron density (occurring 
for the case of y — deg, — 3-db ray, at 77 /X ~ 10 (Fig. 5) extended over an 
arbitrarily large region of some 50 and by neglecting the compensating nega- 
tive gradients. This calculation showed that at most, the curvature could be 
3 X 10 * rad, which is compared with the threshold value at which this effect 
first starts to influence RADVS system performance. This leads one to con- 
clude that the beam will not be significantly bent by the presence of the plume. 
The assumption applied earlier of straight beam rays, is seen to be justified. 

There is one influence of the one-dimensional (plane stratified) plume model 
that requires comment; namely, that while this model precludes bending of the 
beam by the electron density gradients (the rays enter and leave the plume at 
the same angle as long as N e < N,„ the entering and exiting rays will suffer 
lateral displacement. It is straight forward to show that this displacement ( 8 ) is, 
to an adequate level of approximation, represented bv the expression 




Again, considering the worst case ( — 3-db ray, y = 5'- deg) evaluation of this 
expression shows that the displacement will be less than 2 mm, a completely 
negligible amount. 

Consequently, there should be no significant bending or displacement of the 
beam caused by the presence of the plume. 


14 


JPL TECHNICAL REPORT NO. 32-927 


VIL SUMMARY AND CONCLUSIONS 


The Surveyor RAD VS system was analyzed to deter- 
mine if there would be significant interaction between 
the electrons in the vernier thrust chamber exhaust 
plume and the altimeter and doppler beams. Specifically, 
the worst case of interaction, i.e., between doppler beam 
No. 2 and vernier thrust chamber plume No. 1, was 
investigated in detail. It was assumed that the maximum 
specification propellant impurity level (500 ppm) was 
composed solely of sodium. This permitted the thrust 
chamber ionization level to be determined. Frozen flow 
conditions were then assumed to prevail, permitting the 
electron density profile of the exhaust plume to be cal- 
culated. It should be appreciated that these assumptions 
exaggerate the electron density level as the propellant 
impurity concentration is probably on the order of 10 ppm 
and recombination undoubtedly occurs, to a certain ex- 
tent, in the nozzle and the plume (although, for the 
assumed sodium concentration, good correlation was ob- 
tained between theoretical and experimentally deter- 
mined exhaust plane electron density levels). Application 
of the phase integral technique (geometrical optics) per- 
mitted the interaction between the doppler beam and the 
electrons in the plume to be evaluated. The results of this 
evaluation are summarized in the following paragraphs. 

A . Phase Shift 

There is a possibility that two-way (i.e., the transmitted 
beam and the returning beam, neglecting effects attribut- 
able to the characteristics of the surface) nonstatic phase 
shifts of up to 18 deg could occur for various levels of 
thrust chamber gimbal angle and throttling level. The 
threshold value for phase-shift, above which the RADVS- 
system operation would be affected, is — 1 deg. However, 
it should be reemphasized that these calculations were 
performed on the basis of the maximum allowable pro- 
pellant impurity level (500 ppm) being pure sodium. 
Since the phase shift (</>,,) is directly proportional to the 
ionization level c, any reduction in the contamination 
level would be reflected as a decrease in c and a pro- 
portional diminution of <f> p . 

B. Attenuation 

The maximum calculated two-way attenuation is less 
than 10“- db which is well below the allowed 0.1 db. 


C. Reflection 

The operation of the RAD VS system would not be im- 
paired provided that the signal reflected from the plume 
is at least 90 db below the transmitted signal level. 
Under the worst condition of gimbal angle, beam-ray 
position and thrust level, the maximum calculated reflec- 
tion level is 106 db below the transmitted level. 


D. Beam Bending and Displacement 

Under the most exaggerated conditions the calculated 
bending of the beam will be approximately equal to the 
threshold value of 2 mrad, above which the operation of 
the RAD VS system would be influenced. The beam would 
be laterally displaced an amount computed to be less 
than 2 mm, due to its passage through the plume. 


Therefore, it is concluded that, under the most severe 
conditions of geometry and thrust chamber operating 
levels, the phase shift associated with the two-way 
passage of the doppler beam could be of sufficient mag- 
nitude to influence the operation of the RAD VS system. 
The two-way beam attenuation, bending, displacement 
and the back scattering from the plume were found not 
to influence the RADVS system, under the most severe 
operating conditions. 


It should be pointed out that two factors were not 
considered in this analysis; turbulence in the plume and 
the presence of dust particles caused by impingement of 
the plume on the lunar surface. Turbulence could lead to 
a nonsteady increase in the reflectivity of the plume, and 
the presence of dust could lead to an increase in the elec- 
tron density level. Moreover, the analysis is based on a 
one-dimension geometrical-optic treatment of the beam- 
plume interaction, three dimensional effects having been 
neglected. The magnitude of these effects is thought to 
be comparatively small, but they should be evaluated 
as part of any future study. 
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